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Abstract 

It is a well established result that, in classical dynamical systems with sufficient 
time-scale separation, the fast chaotic degrees of freedom are well modeled by (Gaus- 
sian) white noise. In this paper, we present the stochastic dynamical description for 
intermediate time-scale motions with insufficient time-scale separation from the slow 
dynamical system. First, we analyze how the fast deterministic dynamics can be viewed 
as stochastic dynamics under experimental observation by intrinsic errors of measure- 
ment. Then, we present how the stochastic dynamical description should be modified 
if intermediate time-scale motions exist: the time correlation of the noise ^ is modified 
to {S,i't)(,{t')) = C{x,p)6{t — t'), where C{x,p) is a smooth function of the slow coor- 
dinate {x,p), and generally the cumulants of ^ except its average vary as a smooth 
function of the slow coordinates {x,p). The analysis given in this work actually shows 
that, regardless of the sufficiency of time-scale separation, any complex (chaotic and 
ergodic) dynamical system can be well described using Markov process, if we perfectly 
construct the deterministic part of (extended) stochastic dynamics. 



1 Introduction 

In classical many-particle systems with time-scale separation such as Brownian motion, the 
rapid dynamic fluctuations by numbers of fast and chaotic degrees of freedom are well mod- 
eled by Gaussian white noise, and the Langevin equation gives successful descriptions for 
*E-mail address: junchul@kaist.ac.kr 
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such thermodynamic systems [Il|2]. Moreover, recent theoretical works [3], HI O [6] indicate 
that, even for few fast chaotic degrees of freedom, the dynamic fluctuations are well approx- 
imated by suitable stochastic processes with white noise and the stochastic dynamics still 
gives reasonable descriptions. Such results certify that the stochastic process is generally a 
well founded representation for fast chaotic degrees of freedom, if the time-scale separation 
is sufficient appropriately. Actually the stochastic modeling of fast chaotic dynamics based 
on time-scale separation has been applied successfully in diverse areas of science such as 
hydrodynamics [3 [8], climate models [9], and chemical and biological systems [8l[T0]. In 
many cases of real physical systems, however, the corresponding dynamical systems contain 
diverse intermediate time-scale motions and cannot be well decomposed or approximated 
simply by two sub-systems with the sufficient time-scale separation. In this paper, we ad- 
dress the problem of how the intermediate time-scale motions should be represented in the 
point of view of stochastic dynamics. We clarify the role of the intermediate time-scale 
dynamics in the stochastic description of classical Hamiltonian systems and its qualitative 
difference from the fast and slow dynamics which are represented as Gaussian white noise 
and the deterministic part in the Langevin equation, respectively. 
Consider the Hamiltonian 



in which a slow system Xg = {x,p) is interconnected with a fast system Xf = {qi,Pi) = 
(gi,Pi, . . . ,qN,PN) by a potential V{x,qi). Assuming that the time-scale separation is suf- 
ficient, for the motion of p by Hsiow, the potential V{x,qi) introduces additional dynamics, 
which is composed of slow dynamics, such as damping in Brownian motion [llj, and rapid 
fluctuating dynamics. We represent the rapid fluctuating term in the equation of motion 
for p as k{x,p, qi,Pi) and write ^ = —dxHsiow — d^V = h{x,p) + k{x,p, qi-iPi), in which the 
slow dynamics by V is contained in Then, we consider the following deterministic 

system: 



where x = {x,p,qi,pi, . . . ,qN,PN) = {xs,Xf) and i = 1, . . . , N. In the case that we cannot 

^As is well known from the projection operator method [1], in equilibrium thermodynamical system, 
h{x,p) can be extracted by the projection operation, if one can perfectly construct the relevant subspace in 
Hilbert space. 






(2) 
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know or it is unnecessary to know the exact information for the fast variable part in the 
equations of motion ([2]), we may approximate ([2]) to stochastic differential equations, for 
example, 

|=P, f^=hix,p)+^, (3) 

with a suitable time correlation function (^(t)^(t')) ^ Gaussian probability density func- 
tion (PDF) for a random variable ^. In this case, the stochasticity in ([3]) is just a result 
of mathematical approximation or simplification of the original complex system. Such ap- 
proximation can be validated only through the observed empirical data for the system ([2]). 
Here, it should be noted that experimental data always contains unpredictable intrinsic 
errors due to the imperfection of measurement process, which destroy the information of 
fine dynamical structures. The 'unpredictable' nature of error makes such fine dynamics 
change into stochastic dynamics in the resultant experimental data. If we assume all obtain- 
able deterministic information for (1) in our measurement systems is described by the term 
h{x,p), actually the mathematically introduced stochasticity by C, should correspond to the 
inevitable stochasticity in the experimental data originated by the intrinsic measurement 
errors: the stochasticity contained in the experimental data gives the physical identification 
for the mathematically introduced. In what follows, we investigate how the stochasticity 
originated by time measurement errors is related with the fine dynamical structures. 



2 Stochasticity induced by time measurement error 

Let us assume that there are inherent errors in our measurement system for coordinates and 
time. We obtain an erroneous value {x,t) in measurement at an arbitrary instant, while 
the ideal coordinates and ideal time, which is measured by a perfect measurement system 
without error at the same instant, is (x*, t*) = {x, t) + (Ax, At); that is, (Ax, At) is the error 
in our measurement. Under our measurement system, the dynamics of ([2]) is observed as the 
time series x(t), which is deformed from the original dynamics of ([2]), as follows: 

f(t) = x*(t*) + Ax, (4) 

where x*(t*) satisfies the equation of motion ([2]) exactly. In this work, we analyze how the 
dynamics of x(t) changes as the accuracy of our measurement system changes; if we assume 
that, for an arbitrary measurement system, the statistical properties of the error (Ax, At) is 
always independent of time and the error occurrences at different times are probabilistically 
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independent of each other, the time-dependence of x{t) only comes from x*{t*), and the 
coordinate measurement error Ax does not have any role in our analysis — only the time 
measurement error At is important in analyzing the variation of the dynamical behavior 
of x{t) in relation with the change of the accuracy of measurement systemic As can be 
easily checked, all results given in this work essentially hold, regardless of the presence of 
the effect by coordinates measurement error Ax. Thus, to avoid notational complexes and 
to simplify the arguments, we only consider the effect by time measurement error At, which 
does not make any loss of the generality of the analysis. We assume that, except time, we are 
perfectly informed about all values of the observables x, p, qi, and pi in ([2]) at each instant 
from a perfect observer. The situation under our consideration is that there is a missed 
time-parametrization for the perfectly observed x at each instant, i.e., 

x{t) = s*{n, (5) 

but generally 

t ^ t*. (6) 
In the following, we use the notation for x*{t*) as 

= x(t*) (7) 

in order to emphasize that x = x* at an arbitrary instant; we obtain the time series x{t), 
while the perfect observer obtains x{t*), which satisfies the equations of motion ([2]) exactly. 

We can characterize the time measurement errors by giving the PDF ft{At), for which 
the probability to obtain an error in the range {At, At + d{At)) is given as ft{At)d{At). 
Then, we can be confident about the measured data t within an error range as 

t* <t <t* + -, (8) 
2 2' ^ ' 

where we can define the suitable | based on the PDF ft{At). Thus, for two measured times 
ti and t2 such that < ^2 — < we cannot be certain > ti because their error-ranges 
are overlapped, i.e., to ensure that the measured times give the correct time order of t^ and 
t2, the measured times should be separated as 1^2 — ti\ > e. In other words, approximately 
there is a minimum length of time elapse that can be identified by a given time measurement 

^In this respect, time measurement error plays a special role in dynamics. 
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system (one cannot chase correctly the causahty of the events occurring within a time interval 
shorter than the minimum), and this minimum length of time elapse is determined by the 
inherent errors in the measurement system. In translating the deterministic system ([2]) into 
a stochastic system like ([3]) to explain experimental data, an important physical meaning 
is added to the differential dt. The physical meaning of dt in (I3l) is the minimum length of 
time elapse identifiable in our measurement system. We can determine the minimum length 
of time elapse dt to be e providing that the probability to obtain a time value outside of the 
range t* — | < t < t* + | is very small: 



Let us assume that the dynamical system ([2]) starts from an initial value of x. Then, if 
we measure an observable A{x) at t for system ([2]), the measured value of A at t, i.e., A{t) 
is exactly given by the equations of motion ([2]) as A{t*) = A{x{t*)), where t* is a value in 
t-|<t*<t + |: t* randomly has a value in the range with the PDF /t(At). Thus, Ait) 
is given by one of the elements in the following set: 



which is the set of all values obtainable in the measurement of A at t. Exactly which value in 
(3_4(t, e) is given for A in our measurement at t is totally a probabilistic problem determined 
by the following two components: (i) PDF /t(At) and (ii) the set where we assume 

that /f(At) is invariant for time translation and the error occurrences at different times 
are probabilistically independent of each other. In what follows, we assume /t(At) is well 
approximated as /i(At) ^ ^ for — | < At < |, otherwise fti'^t) ~ 0. 

3 Conventional Langevin equation system as a special 
case 

Let A in (ITOl) be the fast coordinates Xf. Then, the set &gj,{t,e) = &(t,e) is the trajectory 
of Xf(t*) during t — |<t*<t + |in the fast phase space F defined hy xj. As the value of e 
increases, i.e., as our measurement system is less fine, the set &{t,e) occupies larger part of 
Q, where Q is the set of all solutions of ([2]) in F. Let us consider the extremal case where 
e is sufficiently large for Xf to wander over almost all areas of Q during t — |<t*<t + |. 



dt = e. 



(9) 




(10) 
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providing that Xf is ergodic and chaotic [T2]. Then, the set &{t,e) is almost the same as Q 
and independent of time, which also guarantees the time independence of the set ©^^(t, e) 
for any fast observable J-'{xf). Thus, if x/ is fast enough to satisfy the extremal case for 
a given e, the statistical properties of the set &{t,e) are almost stationary for time: that 
is, the statistical properties of the measured values of J-'{xf) at t are stationary, and T{xf) 
behaves as a random variable with stationary statistical properties. On the other hand, it 
should be noted that the set &{t,e) is to be slowly dependent on time, because &{t,e) is 
determined by the fast equations of motion from the Hamiltonian 

n' = nfast{q^,P^) + V{x,q{) (11) 

and the value of Ti' has a slow time-dependence through the slow motion of x originated 
from h{x,p) in {XSj. Thus, the extremal case is actually the case where the potential 
V{x, Qi) is regarded as an extremely slow varying function for the deterministic slow motion 
of Xs, i.e., V{x, Qi) is effectively invariant for the slow motion of Xg and the trajectory of x 
is almost confined near the hypersurface Ti' = const (oscillating around the hypersurface 

Also, in the extremal case, by the ergodic property, the time correlation function of any 
fast observable J-'{t*) = J-' {xf{t*)) approximately satisfies 

{HtdHt*2))-0 for \tl-t;\>e, (12) 

where we assume (jF(t*)) = 0. In our measurement system, (concerning the time mea- 
surement errors) the experimentally ascertained quantity as the correlation function of J-'{t) 
is 

{T{ti)T{t2))exp = !\ (^(tl + Ati)^(t2 + At2)) 

f f 

x/i(Ati)/i(At2)rf(Ati)rf(At2), (13) 

where, in the integrand, we use ^(^i 2) = -^(^1,2 + Ati_2). If \ti — t2\ > 2e, we obtain 
{J^iti + Ati)jF(t2 + At2)) ~ from ( fT2l) for any Ati and At2 in the integration ranges, and 

{Hhmh)),., ^ for |ti - t,\ > 2e. (14) 



^Here, we assumed that the oscihation energy is comparatively smaU to the average of H' during dt for 
simphcity in argument, but, without the assumption, (|17p and (|19|) can be derived in similar proof process. 
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Thus, the correlation actually behaves as a 5-function in our measurement system, i.e.. 

Let us analyze the behavior of the rapid fluctuating time series k{x) in ([2]) under our 
measurement systemlf] Expanding k{x) for the slow coordinates Xg, we obtain k{x) = 
JFo(x/) + YlTm.=i'^i-m{,Xf)x^p'^. As prcviously argued, the first term jFo(x/) behaves as a 
white noise with stationary statistical properties in the extremal case. The remaining part 
Afc(x) = k{x) — J-'olxf) also should behave as a white noise in the extremal case, as shown 
in the following argument. Let us denote by (■ ■ ■ )s^(t,(it) the average for all elements of 
©^(t, dt). Assuming that all linear increment of p{t) during dt is given by the deterministic 
part as h{x,p)dt, we request {k{x)) Q^^t^dt) = 0, which means that {Ak{x))eg{t,dt) = 0, because 
{J^o{^f))eit,dt) = const = 0. Thus, we can conclude that the time series Ak{x) is composed of 
the fast fluctuating motion of Xg and the motion of a;/. If we assume that the fast fluctuating 
motion of Xg has the same time-scale as x/, then the motion satisfies a ^-function correlation 
as (fTi|) . and Ak{x) also behaves as a white noise. Thus, we can replace k{x) by a single 
white noise ^, and the equation of motion of p is observed as 

^ = h{x,p)+^, (15a) 

mat')) = 2C6{t-t') (15b) 

in the measurement system with dt{= e), where C is a constant. Also, it is reasonable 
to assume that the PDF for ^, which is determined basically by ft{At) and &{t,e) ~ Q, 
is Gaussian [14j: ^ can be regarded as the fluctuation from an ideal thermal reservoir in 
equilibrium, because the behavior of ^ is independent from Xg and the energy flux between 
the fast and slow system during dt is averagely 0. 



''The time t = and t ^ e respectively corresponds to i = and < ^ in the case where t is continuously 
variable. Since linit^o needs not to be 0, (jl4|) is sufficient for (5-function correlation. 

^The term h{x,p) gives the averaged dynamics for all possible measurement errors; the errors cannot 
deform this deterministic dynamical structure. Basically, we assume that the times series h{x,p) is sufficiently 
slow so that h{t* + e) — h{t*) « h{t*)e, and the error-induced variation in the measured value of h{x,p) is 
comparatively negligible to that of k{x). Thus, we analyze the fluctuating term k{x) only. 
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4 Stochastic dynamical description for insufficiently time- 
scale separated dynamical system 



Now, we generalize the above argument. Let us denote by t the minimum of the time length 
6 for which we can approximate &{t,6) ~ Q for the fast dynamical system. Then, if the 
relative accuracy of our measurement system, denoted by Wloid, to the fast dynamics is not 
fine so that dt{= e) > r, it corresponds to the extremal case. 

Let us observe the dynamical system ([2]) with a finer new measurement system denoted 
by M.new such that dt < r. Then, there is some non-ignorable area in Q which can not 
be reached by Xf during dt, and we cannot simply approximate &{t,dt) ^ Q, i.e., we have 
&{t,dt) C Q. The set &{t,dt) becomes dependent on t, and the statistical properties of the 
set &{t,e) are generally non-stationary: J-'{xf) generally behaves as a random variable with 
no n- stationary statistical properties. Considering &[t,dt) is the trajectory of Xf{t*) during 
t — Y < t* < t + Y, its time-dependence indicates that &{t,dt) is dependent on the value 
of x(t*) = {xs,Xf)(t*) at t* = t. In here, each x^-dependence and x/-dependence has some 
different dynamical meaning. In the following analysis, we will show that the ^/-dependence 
is actually related with the deterministic slow motions which are newly emerged in M.new 
and the x^-dependence is related with intermediate time-scale motions. 

If we expand k{x) in ([2]) for the fast coordinates Xj, we obtain k{x) = J2\a\>o'^a{xs)xJ, 
where a = (ai, . . . , a2N) is a multi- index and we can set Sq{xs) = 0. Then, the linear term 
k{xf) = X]|o|=i '^a(O)^/ generates the most slow dynamics in k{x) — while all xj for |a| > 1 
are detected as white noises in Moid, Qi and generate the most slow time series among 
xj. As our measurement system changes from M.oid to M.new, the variation of the statistical 
properties (the moments or cumulants) of the measured values of k{x) at t mainly comes 
from the most slow time series k{xf). In the following, we analyze how the most slow time 
series k{xf) is viewed in M„ew {dt < r). 

4.1 ^/-dependence of the set <3{t,dt) 

Let us assume that &{t,dt) is only dependent on Xf and almost independent of Xg. As 
previously argued in the extremal case, the Xs-independence of &{t,dt) actually means that 
the motion of x/ is nearly confined on a hypersurface Ti' = const; Q is approximately 
given as the Ti' = const surface. Thus, in this case, the x/-dependence means that dt is 
not long enough for x/ to wander over all areas of the H' = const surface. Consider the 
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averaged dynamics {xf)e{t,dt)- Because Xf cannot wander over all areas of the TC' = const 
surface during dt < r, generally {xf)e(^t,dt) 7^ const for time in Mnew] on the contrary, in 
Mold with dt > T, {xf)e(t,dt) = const = (with a suitable setting of the origin of F), 
because &{t,dt) ^ (5(t, r) ^ Q from the definition of r. In M.oid with dt > r, {xf)Q(^t,dt) 
is stationary as is fixed at the origin of F but, as our measurement system is finer so that 
dt < r, {xf)e[t,dt) is not fixed any more and starts to move slowly with time in f|^ That is, 
the ensemble average of the measured values of xj at each instant varies slowly over time 
in M.new, which is an undetectable deterministic motion in Mo;,^ and viewed as a stochastic 
process in Moid- Considering {k{xf))e{t,dt) = Z]|a|=i '5a(0)(x")s{t,dt) also slowly varies by 
the motion of {xf)Q(^t,dt), eventually, this newly emerged slow dynamics induces the observer 
using M„e^ to change the deterministic term h{x,p) in ([3]) to h{x,p) + 6h{x,p) by adding 
the more detailed information of the dynamics. 

On the other hand, the motion of the fast fluctuating part 

yf = ^f- {^f)eit,dt) (16) 

is confined on the energy surface H' = const in the yf phase space (note that {xf)e(t4t) is 
a function of x/), and the set &yj{t, dt) contains all points of the Ti' = const surface, under 
the assumption that the motion of yf is still ergodic and chaotic: if yj cannot wander over 
all areas of the energy surface during dt, generally {yf)6{t,dt) 7^ const, which contradicts with 
the definition of yj. Thus, the motion of yj becomes the extremal case again in Mnew with 
dt < T. 

There is just a replacement of the deterministic part h{x,p) by h{x,p) + 5h(x,p) and 
the fluctuating noise part is identical to the case of the Langevin equation or the extremal 
case, i.e., white noise with stationary statistical properties; the x/-dependence of the set 
&{t, dt) changes the time series k{xf) from white noises in Moid to colored noises in Mnew as 

^{^f)e(t,dt) is a function Xf from the ^/-dependence of &{t,dt). 

^Strictly speaking, Sh{x,p) is over-simplified expression for the newly emerged slow dynamics. If we 
define a new slow variable x' = (fc(a?/))s(t,tit)i then ^ = h{x,p) + 2:' -f X]|q|=i SaiS^)y^- Generally, x' may 
be independent of x and p in the sense that, for some ti ^ t2, {x,p){ti) = {x,p){t2) but x'{ti) 7^ x'{t2). 
We have to contain the new additional equation of motion in the deterministic part as (in the most general 
expression) ^ = hnew{x,p,x' ,t). But, x'{t) is very slow varying time series, and there can be one-to-one 
correspondence between {x,p,x') and t for very long time interval. Thus, we can write the deterministic 
part of the equations of motion as ^ — p, ^ — hnew{x,p, x'), and ^ = h{x,p) + x' for very long time. 
Under M„e^, generally the number of the slow variables may be increased, e.g., Xg = {x,p, x')] thus, in our 
notation, we can regard {x,p) under M„eu, as a more expanded slow variable set than {x,p) under Mo;^. If 
we intend to write the equations of motion only using x and p except x', generally the equations of motion 
have to contain the memory term for {x,p) instead of x' . 
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the combination of the deterministic motion 6h{x,p) and white noises. 
4.2 Xg-dependence of the set G{t,dt) 

Next, we analyze the af^-dependence of the set G{t,dt). Let us assume that &{t,dt) is only 
dependent on Xg and almost independent of Xf] the V{x, qi) is not effectively invariant for the 
slow motion of Xg any more. Actually, the x/-independence means {xf)e(^t4t) = const = 
for time, as shown in the following argument. If we fix the value of Xg arbitrarily, the set 
&{t, dt) is invariant for any value of on the hypersurface Ti' = const corresponding to the 
fixed value of Xg (note that Ti is conserved), which means that the set &{t,dt) contains all 
points Xf on the hypersurface. Also, since Ti' varies slowly with the slow motion of Xg, the 
elements of &{t,dt) are confined near the hypersurface. Thus, actually the average for the 
set &(t, dt) can be approximated as the average for all points on the H' = const surface and 
{xf)6{t,dt) = 0, which holds for arbitrary values of Xg. 

The Xg-dependence induces the change of the set &{t,dt) while keeping {xf)e(^t,dt) = 
invariant for time, which means the following two results: 

(i) the variation of Xg induces the variation of the other statistical properties of &{t,dt), 
such as the averages of qf, qiPj, qiqjPk, • • • , which are stationary in Moid, i-e., generally 
the statistical properties {xj)e(^t,dt) for |a| > 2 are non-stationary and vary as a smooth 
function of Xg in M„e^„I^ and 

(ii) the Xs-dependence does not make any change in the deterministic part h{x,p) in ([3]). 

Except {k{xf))e{t,dt), generally {k°'{xf))eit,dt) for 1^1 > 2 vary as a function of Xg in Mnew, i-e., 
the statistical properties of the measured values of k{xf) at t, except the ensemble average 
of k{xf), vary slowly function of Xg in M„ei«- 

Also, the time correlation function of k{t*) is given as {k{tl)k{tl)) ^ for |t* — t^l > dt, 
because dt is sufficient for Xf to wander over all areas of an energy surface — the &{t,dt) 
contains all xj on the energy surface. Thus, considering {k'^{xf))Q(t^dt) should be a smooth 
function of Xg and using (|T3l) and f|T^ . the motions inducing the Xg-dependence give the 
following time correlation for k{t) in M„ew with dt < r: 

(e(tl)e(t2))e.p = C(x,p)5(ti-t2), (17) 

^The smoothness of {xj)e(t.dt) ^ function of Xs comes from the fact that hmt^^t^ &{ti, dt) ~ 6(^2, dt) 
and the trajectory of Xf in F is smooth. 
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where the notation k{xf) is replaced by ^ and C{x,p) is a smooth function of Xg. Therefore, 
under the x^-dependence, k{xf) still behaves as a white noise as satisfies (fTTl) . but, differently 
from the conventional white noise with stationary statistical properties, its moments or 
cumulants except its average are non- stationary for time as a smooth function of the slow 
coordinates Xs- 



4.3 Extension in Stochastic dynamical formalism 

As is clear from the above analysis, the motions of x/ generating the a;<i-dependence have 
a faster time-scale than the motions generating the x/-dependence which are observed as 
the deterministic slow dynamics 6h{x,p): under the ^/-dependence, xj is not fast enough 
to wander over all areas of an energy surface during dt but, under the x^-dependence, xj is 
fast enough to do so. Also, the motions of Xf generating the x^-dependence have a slower 
time-scale than the motions of yf = xj — {xf)Q{t,dt) in the case of x/-dependence, i.e., the 
fast system in the extremal case (the fast system described as white noise in the Langevin 
equation): under the Xg-dependence, the motions of x/ cannot be approximated as confined 
near an energy surface, and, since the elements of &{t,dt) are confined near an energy 
surface, x/ cannot wander over the entire areas of Q during dt. Thus, actually the motions 
related with the Xg-dependence are intermediate time-scale motions. 

Consequently, the previously obtained results for the case of x^-dependence show that, 
if there exist intermediate time-scale motions, generally the fast dynamics, which cannot be 
written in the deterministic part of stochastic dynamics, should be described as the white 
noise ^ of which the cumulant 



dk"- ^ m\ 

m=0 



fc=0 



for n > 2 is a smooth function of the slow coordinates x^E Thus, we have to treat the 
cumulants c„(x,p) as dynamical variables in stochastic dynamics and have the following 
stochastic dynamics for the coordinates (x,p, C2, C3, . . . ) with the setting c\ = 0: 

dx dp , . . ^ /-.^ N 

-^=P, -^ = Kx,P)+i, (19a) 



'The PDF for ^ is a function of Xg, and generally the functional form of the PDF varies as Xs varies. 
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^ = + ^± ^hixpn (i9h) 

dt dx dt^ dp dt-^""^^'^'^'' ^^^> 

where n = 2, 3, . . . , and 

mm) = c2{x,p¥{t-t') (20) 

by the result (fT7|) : the cumulants c„ and (!20|) completely determine any time correlations 
(^"K'^i) ■ ■ where rii, . . . ,ni are positive integers, and thus the formulation f|T9l) to- 

gether with ( l20l) gives a complete stochastic dynamical description. In deriving ( |T9l) . sim- 
ilarly as in the extremal case, k{x) can be replaced by a single noise C,'- assuming that all 
linear increment of p{t) during dt is given by the deterministic term as h{x,p)dt, we have 
{Ak{x)) Qg{t,dt) = 0, where Ak{x) = k{x) — k{xf), and, also if we assume that the fast 
fluctuation of Xg has the same time-scale as x/, Ak{x) should satisfy at least (|T7j) . How- 
ever, generally the time-scale of xj becomes faster as |a| increases, and, for some |a| = M, 
the term J2\a\>M '^a{xs)x'j behaves as a white noise. Thus, more precisely, we can write as 
^ = h{x,p)+^+ri, where rj is the conventional Gaussian white noise and {^{t)r){t')) oc 6{t—t'). 
While the trivial case (c2 = const ^ and c„ = for n > 3) of formulation ( fT9l) gives the 
conventional Langevin equation, in the most simple nontrivial all c„ = except 

C2 7^ const, (as shown in the subsequent argument) the formulation can be reduced to the 
conventional multiplicative noise method. 

The results obtained so far can be summarized as follows. Intermediate time-scale mo- 
tions cannot be described by deterministic trajectories in phase space, because they satisfy 
a (5-function time correlation as ( |T7I) . and also they cannot be described simply by random 
forces acting on the trajectories, because they have deterministic properties represented by 
the cumulants Cn{x,p) for n > 2. Intermediate time-scale motions should be described by the 
trajectories in the cumulant space defined by the coordinates (c2, C3, . . . ). This indicates that 
the simple trajectory-based description in phase space, such as the Langevin or generalized 
Langevin equation, cannot describe systematically the dynamical effect from diverse inter- 
mediate time-scale motions. Also, since intermediate time-scale motions satisfy a (5-function 
time correlation, it means that, regardless of the sufficiency of time-scale separation, any 
complex ( chaotic and ergodic ) dynamical system can be well described using Markov process, 
if we perfectly construct the deterministic part in stochastic dynamics. 

In formulation (|T9l) . ^ and {x,p) exchange their influence with each other. Thus, the 
random noise ^ cannot be considered as the thermal fluctuation from an ideal thermal 
reservoir in that the statistical properties of ^ are not independent from the system described 
by {x,p) — the thermal reservoir by the fast chaotic dynamics of Xf is treated as a finite 
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thermal system not having infinite capacity. 

In the most simple case of formulation ( |T9l) as all = except C2 = 2C{x,p), we have 
(^(^)^(^')) = 2C{x,p)5{t - t') from ([20]), and the PDF w(0 for ^ is given as a Gaussian: 
'w{^) oc e^^^/^''*^'^''^^ If we define C, = ^/C{x,p)^ , the variance of the random variable ( is 
constant, and, using (, we can reformulate the most simple case of (fT9l) in a conventional 
multiplicative noise system, as follows: 

| = M.,p)-i^ + C(.,rtk, (21a) 

(CWC(i')) = 2i(«-«'). (21b) 

where the equation of motion for x is the same as in (fT9l) . and the PDF for ( is Gaussian: 
(C) oc e-f'/^ As can be easily checked, these two systems give a same Fokker-Planck 
equation; in this respect, the multiplicative noise method is equivalent to the most simple 
case of formulation ffT9l). 



5 Additional remarks 

Finally, we point out that the arguments for intrinsic time measurement errors symmetri- 
cally hold for space measurement errors in static problems, e.g., dxip{x, t) = hs{x, t) + hf{x, t) 
where hs and hf are slow and fast varying parts for x, respectively. There is the experimen- 
tally identifiable minimum distance dx corresponding to dt (in dynamics), and we obtain the 
stochastic statics, dxip{x,t) = hs{x,t) + ^(x) for a random variable ^(x). The stochasticity 
originated by space measurement errors refiects the fine structure of statics on the statistical 
properties of ^{x). 
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